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Abstract 

A theory is presented for the modification of bandgaps in atomically thin boron nitride (BN) by 
attractive interactions mediated through phonons in a polarizable substrate, or in the BN plane. 
Gap equations are solved, and gap enhancements are found to range up to 70% for 
dimensionless electron-phonon coupling A = 1, indicating that a proportion of the measured BN 
bandgap may have a phonon origin. 

Introduction 

The need for bandgaps in graphene on eV scales has led to a number of proposals, such as 
the use of bilayer graphene [I] creation of nanoribbons [2] and manipulation through 
substrates [3j[4]. Recently it has become possible to manipulate atomically thin layers of 
boron nitride (BN) and other materials with structure similar to graphene [5). This may 
lead to a complimentary method of manipulating bandgaps to make digital transistors. 
In low dimensional materials, strong effective electron-electron interactions can be induced 
via interaction between electrons confined to a plane and phonons in a polarizable 
neighboring layer |6]. Theory has shown that similar interactions account for the transport 
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properties of graphene on polarizable substrates [7], and that sandwiching graphene 
between polarisable superstrates and gap opening substrates can cause gap 
enhancement pi. This paper examines similar gap changes in atomically thin BN due to 
interactions mediated through substrates. 



Model 

Atomically thick hexagonal BN (h-BN) has similar chemistry to graphene: Bonding occurs 
through sp2 hybridization, and electrons with energies close to the chemical potential are in 
unhybridized 7r orbitals [9]. A key difference is that the electronic charge is not completely 
screened by the sp2 hybridization, shifting tc orbitals by A n = +A on N sites and —A on B 
sites. This shift is the dominant cause of a gap of order 2 A. Tight binding fits to ab-initio 
simulations of monolayer BN have established that the hopping, t = 2.33eV [10] , with an 
estimate of A = 1.96eV= 0.84t. Experiments indicate larger gaps: bulk h-BN has 



5.971eV 111], and monolayer h-BN has a gap of 5.56eV 12 corresponding to 
A = 2.78eV= 1.20t. There is significant variation in phonon energies, hQ q , in h-BN 13 
LA phonons energies range up to around 140meV at the M point, and TA phonons to 
around HOmeV at the K point. Optical phonon energies range between 160meV and 
200meV. Coupling, f n (m), between electrons and phonons in either a polarisable 
substrate, or the BN monolayer is possible and the corresponding Hamiltonian is, 

H = -t ^ ( a na c ria + c ] nl(7 a na )- ^ fn(m)n nrT ^ m + ^hn(N Tn + l/2)+^2 ^ n n na . (1) 

(n,n')a nma m na 



The Hamiltonian terms are shown schematically in Fig. [I] (left). a) n(J creates electrons of 
spin a on B sites and c^, CT on N sites. Vectors n are to atoms in the monolayer, and m to 
atoms in the substrate. N m is the number operator for phonons. The Hamiltonian is also 
approximately valid for interactions in the plane, and Fig. [I] (right) shows the forces on 
ions from an increase in electron density at a B site. The largest forces are on the 
near-neighbor sites, so the effective interaction is mainly site diagonal (electrons on A sites 
self-interact through phonons on B sites and vice versa). The diagram indicates that the 
strongest interaction is between electrons and optical phonon modes. 
For simplicity, I use the Holstein electron-phonon interaction, f n (m) oc 5 n<n >, which 
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qualitatively captures the physics. There may be quantitative changes to the results for 
longer range Frohlich interactions and from modulation of the electron-phonon interaction 
due to incommensurability of the substrate, which was estimated at around ±8% of the 
average value [8~1. 



Gap equations and results 

Low order perturbation theory is applicable for low phonon frequency and weak coupling. I 
derive a set of gap equations by symmetrizing the self energy, 



v , ihu n (l - Z£) + A^ 



The local approximation used here is a good starting point here because the modulated 
potential A is large, and electrons are well localised. Off diagonal terms do not feature in 
the lowest order perturbation theory for the Holstein model since the interaction is site 
diagonal. Z n is the quasi-particle weight and A n is the gap function. For bosonic 
quantities, hu s = liskBTs and for fermions hu n = l-KksTin + 1/2). T is the temperature 
and n and s are integers. 

The full Green function can be established using Dyson's equation 
G~ 1 (k,iu n ) = GQ 1 (fc,ia; n ) — X(z"u; n ), leading to, 



>k 
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Substituting the expression for the Green function into the lowest order contribution to the 
self energy, 

E u (k,iuj n ) = -ksTtXiY^ [ P L G ii (k-q 1 iu n „ s )(d% i \q,u s )-2d { j i \0,0)). (3) 

^ J vbz 

Here, the phonon propagator, d^\q,oj s ) = 5ijVL 2 /(VL 2 + oj 2 ), hence there are no off diagonal 
elements of the lowest order self energy. The use of a single averaged Q and A here is 
consistent with a mean-field approximation. At half-filling, it is reasonable to assume that 
Aa — within the scope of the model, so that A A = A B and Z A = Z B . This leads to the 
gap equations, 

= E J d %^ L , ZL , + A^ + S- (4) 



t\k B T y-y f D(e)u n _ s Z n _ s d (iuj s ) 



hu n ^fj h A u 2 n _ s Z 2 n _ s + A' n 2 _ s + e 
where the full gap is A^ = A n + A. The density of states for a tight binding hexagonal 



lattice in the absence of a gap, D(e), has the form given in Ref. 14 . The equations may be 
solved self consistently by performing a truncated sum on Matsubara frequencies. 
Gap and quasi-particle weight functions only have a weak Matsubara frequency 
dependence (< 0.3% for A = 1, k^T = HQ = O.Olt). The local gap enhancement factors 
A'/A, are shown in Fig. [2] (left) for various A, showing a modest increase of ~ 70% for 
A = 1. The enhancement factor increases slightly with decreasing A but is essentially 
unchanged by modifications to phonon frequency and temperature for the parameter values 
used here. I also calculate the temperature dependence of the gap, Fig. [2] (right). For very 
large temperatures, where k B T approaches A, there is a drop in the gap size. For 
T < 0.3t ~8000 K, this levels off, and the gap becomes relatively constant. 

Summary and conclusions 

I have presented a theory for the modification of BN band-gaps by interaction with 
phonons. It is of interest to make comparison between the bandgaps of bulk h-BN, 
nanotubes, monolayer h-BN and the theory presented here. Measured bandgaps of bulk 



h-BN are of between 5.8eV |15 and 5.971eV 11 , indicating that interaction between layers 



increases the bandgap, consistent with the theory here. The bulk gap is also higher than 



that for nanotubes (5eV) [16]. On the other hand, Song et al. 12 claim that the gap is 
reduced as BN thickness increases. The above discussion is presented with the caveat that 
the theory requires that hopping between the substrate and the BN monolayer is small. 
Interlayer hopping will affect the bandwidth and bandgap, and direct Coulomb interaction 
with strongly ionic substrates could also affect the band structure if the charge density at 
the surface of the substrate varies dramatically. 

It is also of interest to estimate the magnitude of the bandgap modification due to 
electron-phonon interaction in isolated monolayers of BN. Ab-initio calculations have 
attempted to quantify the magnitude of the interaction between electrons and acoustic 
phonons for small momentum excitations [17] . Extrapolating the interaction, and taking a 
mean- field average (assuming mean momentum magnitude of 4n/9a), the electron-phonon 
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coupling can be estimated as A = (Air /9) 2 E 2 /2a 2 tMQ 2 } taking E\ = 3.66eV from Ref. 17 



M « 12.5 amu, a = 2.5A. The mean energy of longitudinal acoustic phonons lies in the 
range 50 — 75meV, giving a range of A = 0.05 — 0.12, so the contribution of phonons to the 
bandgap is estimated as 3 — 7%. I would expect BN to have stronger interaction with 
optical phonons, since the pattern of distortions around an electronic defect is consistent 
with optical modes (see Fig. 1). 

The BN gap is too wide for digital applications. Recently, it has become possible to 
manufacture silicene, an atomically thick layer of silicon with similar properties to 



graphene 18 , so it may be possible to make GaAs or A1P analogues to BN. Smaller gaps 
could be available from those materials, which might be used to create tunable bandgaps 
for atomically thick transistors. 
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Figure 1: (left) BN-substrate system annotated with interactions. Electron-phonon inter- 
actions between the BN layer and substrate are poorly screened, and large interactions of 
strength f n (m) are possible. Ions in the substrate oscillate with frequency Q. N sites have 
energy +A and B sites —A, opening a gap. The attractive phonon mediated electronic 
interaction / binds electrons onto the same site, effectively enhancing the gap. (right) In- 
teractions in a monolayer of BN. 




Figure 2: (left) Modification of the BN bandgap. The gap enhancement depends mainly on 
A, is weakly dependent on A and shows almost no change with Q. Calculations are made for 
A = t corresponding to a BN gap of 2 A = 4.66eV, A = 1.20t (2A = 5.6eV), and A = 0.84t 
(2 A = 3.92eV, the tight binding fit from Ref. [IS]), t = 2.33eV, HQ = OMt = 23meV, 
m = Omt = 70meV, HQ = 0.06t = 140meV, TiVL = 0.09t = 210meV, covering the full range 



of phonon frequencies in Ref. 13 . /c#T = O.Olt (T = 268K) and A < 1 (right) Variation of 
the gap with temperature, hfl = 0.05t = 117meV and A = 0.2. There is a weak temperature 
dependence, due to the large A, consistent with the measurements in Ref. [15], with the gap 
starting to close only for extremely high temperatures T > 8000i^, presumably above the 
melting point of the material. 
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